A dual iterative substructuring method with a penalty term

نویسندگان

  • Chang-Ock Lee
  • Eun-Hee Park
چکیده

An iterative substructuring method with Lagrange multipliers is considered for the second order elliptic problem, which is a variant of the FETI-DP method. The standard FETI-DP formulation is associated with the saddle-point problem which is induced from the minimization problem with a constraint for imposing the continuity across the interface. Starting from the slightly changed saddle-point problem by addition of a penalty term with a positive penalization parameter η, we propose a dual substructuring method which is implemented iteratively by the conjugate gradient method. Performance of such a dual iterative substructuring method is directly connected with the condition number of a relevant dual system. For η = 0, the proposed method is reduced to the FETI-DP method. For the preconditioned FETI-DP with the optimal Dirichlet preconditioner, it is well-known that the condition number is bounded by a polylogarithmic factor: (1+ log(H/h)) in two dimensions and (H/h)(1+ log(H/h)) in three dimensions. To the contrary, in spite of the absence of any preconditioners, it is shown that the proposed method is numerically scalable in the sense that for a large value of η, the condition number of the resultant dual problem is bounded by a constant independent of both the subdomain size H and the mesh size h. We deal with a computational issue and present numerical results. Furthermore, we extend the proposed method to the three dimensional problem.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

A Domain Decomposition Method Based on Augmented Lagrangian with a Penalty Term in Three Dimensions

In our earlier work [4], a dual iterative substructuring method for two dimensional problems was proposed, which is a variant of the FETI-DP method. The proposed method imposes continuity on the interface by not only the pointwise matching condition but also uses a penalty term which measures the jump across the interface. For a large penalization parameter, it was proven that the condition num...

متن کامل

An Iterative Substructuring Method for the Discretized Stokes Equations by a Stabilized Finite Element Method

Abstract. A simple algorithm of iterative substructuring method as the same way of elasticity problem is proposed for a discretized Stokes equation by P1/P1 element and penalty stabilization technique. Owing to the stability term, solvabilities of local Dirichlet problem, of local Neumann problem for preconditioner, and of the coarse space problem are ensured. Conjugate gradient method with pre...

متن کامل

Iterative Substructuring Methods for Spectral Element Discretizations of Elliptic Systems. II: Mixed Methods for Linear Elasticity and Stokes Flow

Iterative substructuring methods are introduced and analyzed for saddle point problems with a penalty term. Two examples of saddle point problems are considered: the mixed formulation of the linear elasticity system and the generalized Stokes system in three dimensions. These problems are discretized with spectral element methods. The resulting stiiness matrices are symmetric and indeenite. The...

متن کامل

Preconditioning of Saddle Point Systems by Substructuring and a Penalty Approach

The focus of this paper is a penalty-based strategy for preconditioning elliptic saddle point systems. As the starting point, we consider the regularization approach of Axelsson in which a related linear system, differing only in the (2,2) block of the coefficient matrix, is introduced. By choosing this block to be negative definite, the dual unknowns of the related system can be eliminated res...

متن کامل

Primal and Dual Interface Concentrated Iterative Substructuring Methods

This paper is devoted to the fast solution of interface concentrated finite element equations. The interface concentrated finite element schemes are constructed on the basis of a non-overlapping domain decomposition where a conforming boundary concentrated finite element approximation is used in every subdomain. Similar to data-sparse boundary element domain decomposition methods the total numb...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Computers & Mathematics with Applications

دوره 64  شماره 

صفحات  -

تاریخ انتشار 2009